Spin entanglement of two spin-^ particles in a classical gravitational field 
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Abstract 

The effect of a classical gravitational field on the spin entanglement of a system of two spin-^ particles 
moving in the curved spacetime is discussed. The system is described by a two-particle Gaussian wave 

cN : 

packet represented in the momentum space and both the acceleration of the system and the curvature of 
O ■ 

. the spacetime cause to produce a Wigner rotation acting on the wave packet as it moves along a path in 

■ the curved spacetime. By calculating the reduced density operator at a final point, we focus on the spin 

entanglement of the system. In a spherically symmetric and static gravitational field, for example a charged 



black hole, there can be particular paths on which the Wigner rotation is trivial and so the initial reduced 
density matrix remains intact. This causes the spin entanglement to be invariant during the motion. The 



spin entanglement descends to zero by increasing the angular velocity of the mean centroid of the system as 
^ ' well increasing the proper time during which the centroid moves on its circular path around the center. 
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I. INTRODUCTION 



The field of quantum information has made rapid progress in recent years. It is important to 
study all those processes that might have an effect on quantum entanglement. Special relativistic 
effects on quantum entanglement and quantum information is investigated by many authors 
^]. It is shown that though Lorentz transformations can change the entanglement of the spins of 
massive particles, the total entanglement is invariant for inertial observers. Some authors studied 
the entanglement between two modes of a free Dirac field as seen by two relatively accelerated 
observers ]{], Q] • This entanglement is degraded by the Unruh effect and asymptotically reaches a 
non-vanishing minimum value in the infinite acceleration limit. Furthermore, using the formalism 
of quantum fields in curved spacetimes, extension of the quantum entanglement and teleportation 
to the gravitational field in four-dimensional spacetimes, as well in higher dimensional black hole 
spacetimes is discussed [ll . 



121 ] . There exists another approach for discussing the effect of gravita- 
tional field on quantum information in which, using the concept of local inertial frames, the special 
relativistic considerations is extended to general relativity. Using this approach and considering a 
wave packet for the system, Terashima and Ueda discussed a mechanism of spin rotation caused by 
spacetime curvature for spin-^ particles moving in a gravitational field [ijj. It is shown that this 
effect gives rise to a spin entropy production that is unique to general relativity. This means that 
even if the spin state of the particle is pure at one spacetime point, it becomes mixed at another 
spacetime point. In the same manner, the gravitational spin entropy production for particles with 
arbitrary spin moving in a curved spacetime is discussed [bj ]. 

In the present article, regarding the latter approach, we discuss the effect of a gravitational field 
on the spin entanglement of a system consisting of two spin-^ particles moving in the gravitational 
field. For theoretical reasons, it may be interesting to quantify the entanglement for such a sys- 
tem. We describe the system by a two-particle wave packet represented in the momentum space 
with a definite width around the momentum of a mean centroid. We consider the mean centroid 
for describing the motion of the system in the gravitational field which is a classic field in our 
approach. Thus we require to assume that the curvature of the gravitational field does not change 
drastically within the spacetime scale of the wave packet. Using the concept of local inertial frames 
it turns out that both the acceleration of the centroid and the curvature of the gravitational field 
cause to produce a Wigner rotation that affects the wave packet. We confine the argument to a 
static spherically symmetric gravitational field such as a charged black hole which because of the 
mathematical form of its metric, allows us to present a detailed argument. At a given initial point 
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in the curved spacetime, the momenta of the particles are taken to be separable, while the spin 
part is chosen to be entangled as the Bell states. By calculating the reduced density operator, we 
focus on the spin entanglement of the system as moves on a circular path around the center. In 
the present problem, one can more closely investigate how a central massive charged body affects 
a quantum communication that takes place around it. 



II. WIGNER ROTATION FOR TWO-PARTICLE STATES MOVING IN A GRAVITA- 
TIONAL FIELD 

In a curved spacetime the curvature causes to break the global rotational symmetry. Therefore, 
the spin of a particle in general relativity can be defined only locally by switching to an inertial 
frame at each point and then invoking the rotational symmetry of the local inertial frame. A local 
inertial frame at each point of a curved spacetime is introduced through a tetrad e a ^(x) defined by 

e a ^(x)e b u {x)g^ v {x) = r] ab (1) 

where g^ v {x) is the metric that describes the gravitational field and r? a & = diag(— 1, 1, 1, 1) is the 
Minkowski metric jlfj. Greek indices run over the four general coordinate labels, while Latin indices 
run over the four inertial coordinate labels 0,1,2,3. Now, a spin-^ particle in a curved spacetime 
is defined as a particle whose states furnish the spin-^ representation of the local Lorentz group. 



A one-particle momentum eigenstate is specified by \p, a) where p = {\/\p\ 2 + w- 2 c 2 , p) is the four- 
momentum eigenvalue as measured in the local inertial frame and a denotes the z-component of the 
spin. Then, the wave packet in the momentum representation for a system of two non-interacting 
particles, observed in a local frame located at an initial point x^ , can be written as 

\^ (i) ) = Yl J J d 3 p 1 d 3 p 2 ^ aia2 (p 1 ,p 2 )\p 1 ,a 1 ;p2,a 2 ) (2) 

(T1CT2 

where p\ and p 2 are the four-momentum of particles, and ^ aiU2 (p l ,p 2 ) are wave functions deter- 
mining momentum and spin distribution, normalized as 



1 1 d 3 Pid 3 p 2 i*wPi,p 2 )i 2 = i ( 3 ) 

01(72 



provided that (jp' x , a[;p 2 , cr 2 |pi, ai;p 2 , a 2 ) = S 3 (p[ - Pi)<5 3 (p' 2 - P 2 ) S a ' i<T1 S^. We can use 
ty aicr2 (pi, p 2 ) for specifying the momentum entanglement, the spin entanglement or even the en- 
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tanglement between spins and momenta. 

It is required to assume that the spacetime curvature does not change drastically within the 
spacetime scale of the wave packet. Otherwise, we encounter some problems in a classical approach 
for the gravitational field. Suppose that the two-particle wave packet ([2]) has a mean centroid that 
in a semiclassical approach one can regard it to describe the motion of the center of mass of the 
two-particle system. The momentum of each particle is assumed to be distributed properly around 
the momentum of the centroid. Assume that the centroid moves along a specified path x^(t) in 
the curved spacetime. So, the four-momentum of the centroid as measured in a local inertial frame 
will be as q a {x) = e a ^(x)(mdx fM /dT). Suppose that at r, the centroid locates at an initial point 
Xi = x^(Ti). Then, the centroid moves along the path x^(t) and at a proper time r/ it reaches to a 
final point Xf = x^(rf). Consider two adjacent points on the path such that dr be the infinitesimal 
proper time between them and at each point a local inertial frame is located. It is shown [2] that 
these frames are related by a local Lorentz transformation A a b (x) = 5 a b + X a b (x)dT, where 

X a b (x) = \ [a a {x)q b {x) - q a (x)a b (x)} + u» {x)[e b v {x)V ^e\{x)] (4) 

mc z 

where u^(x) is the four-velocity of the centroid and a a {x) = e a ^{x)[u u {x)V u u^ {x)] is its four- 
acceleration as measured in the local frame. A classical force is needed to produce this acceleration, 
so for geodesic motions that a a (x) = 0, no force is needed. The first part of @ is accelerated 
related, existing even in special relativity, and the second part is curvature related, that arises 
from the change in the local inertial frame along the path. One can construct a glob al Lorentz 



transformation A(xf,Xi) by adding the successive local Lorentz transformations as |16l | 



A(x/, x^ = Texp 



/ 



A(x(r))dr 



(5) 



where T is the time-ordering operator and X(x) is a matrix whose elements are given by (jU). 

Now imagine that when the centroid is at the initial point Xi, the corresponding wave packet 
as viewed from a local inertial frame located at Xi is given by ([2]). When the centroid reaches to 
the final point xt, the corresponding wave packet, as viewed from a local inertial frame located at 
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that point, can be written as 



|*W> = U(A 1 (x f ,x i ))(g ) U(A 2 (x f ,x i ))\^) (6) 



^ f 3 / (Aipi)° /■ ,3 / (A2P2)° , T , , , 



Pi J V ^2 
D a , iai (W(A 1 ,p 1 ))D a , 2a2 (W(A 2 , P2 ))\A lPl ,a[-,A 2P2 ,cj' 2 } , 



Cl<T2 

X 



where U(A(xf,Xi)) is a unitary operator, W(A,p) is the Wigner rotation operator corresponding 
to A(xf,Xi) and D a / a denotes the two-dimensional representation of the Wigner rotation operator 
[it! ]. The Wigner rotation operator can be expressed explicitly as [u]] 



'* w(x(T))dr\ , (7) 



W(A(x f ,Xi),p) = Texp 

where w is a matrix whose elements are given by 

w*i (x) = A*j (x) + , ( 8 ) 

p u (x) + mc 

where i and fc run over the three spatial inertial frame labels (1,2,3). Note that wa = 0. 

In the following, as the background gravitational field, we consider a static spherically symmetric 
spacetime described with the metric 

ds * = _ c 2 e 2A(r) dt 2 + e 2B { r) ^2 + r 2 ^2 + ^2 fl # 2) (g) 

The asymptotic flatness assumption imposes the following conditions on A(r) and B(r), 

lim Air) = lim S(r) = 0. (10) 

r— >oo r— >oo 

We require to introduce a local inertial frame at each point of the spacetime. Hence we employ 
the tetrad 

eo^-e"^), ef = e- B <r), e 2 e = -, = (11) 

c r r sin a 

with all the other components being zero. 

Generally, in a given gravitational field, the mathematical form of the local Lorentz transforma- 
tion and the corresponding Wigner rotation depends on the path along which the centroid (wave 
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packet) moves. For the metric ([9]), regarding its spherical symmetry, we choose conveniently the 
equatorial plane 9 = ^ as the plane of motion and suppose that the centroid is moving with a 
constant speed v on a specified circle of radius r around the center. This motion is generally a 
non-geodesic motion and so a central force is needed to maintain the system on the orbit. Then, 
the components of the four-momentum of the centroid in the local inertial frame are 



q° = 7171c 



I 1 = q 2 = 0. 



q = jmv. 



(12) 



In this case, the acceleration has only one non-zero component as 

1 7 2 -1 



i 1 (r) = c 2 7 2 e- B W (A'(r) 



r 7 Z 



and then the matrix A has only four non-zero elements as 

A 1 o(r) = A° 1 (r)=c 7 (7 2 -l)e^» 
A x 3 (r) = -A 3 i(r) = -^ve~ B{r) 



A'{r) 
A'(r) 



(13) 



(14) 
(15) 



These components as substituted in ([8]), lead to the following non-zero time- independent elements 
for the local Wigner rotation matrix w 



w 1 3 



-w 3 i 



-B(r) 



-ce 



mc V \ mc 



A\r) - - 
r 



(16) 



+ 1 



IP 



mcj vP 2 + m 2 c 2 + mc 



Consequently, for circular motions in a general spherically symmetric spacetime, the Wigner rota- 
tion operator ([7j) reduces to 



W(A(t)) 



cos sin 



10 

y — sin cos J 



(17) 



where = w 1 2,t with r = r/ — t s . This operator shows a rotation about the 2-axis and has a 
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2-dimensional representation as 



D(Q) = e 



-ij 2 e _ 



cos y — sin y 



(18) 




COS y 



where J 2 is the 2-component of the angular momentum operator. 

One may argue that if A'(r) = 1/r identically, © will vanish and there is no Wigner rotation. 
However, this condition leads to a logarithmic r-dependence for A(r), which contradicts the asymp- 
totic flatness assumption we made for the metric Q. Of course, the condition can be fulfilled at 
some distinct radii. 

III. REDUCED DENSITY MATRIX 

We intend to calculate the entanglement between the spins of particles described by the final 
state ([6]). To do this we refer to the corresponding reduced density operator obtained by taking 
the trace over the momentum of the density operator pf = \"^f^)(^>^'\, that is 



To follow the argument we need to choose an explicit form for ^f aia2 (p 1 , p 2 ). We assume no 
entanglement between spins and momenta, as well as no entanglement between momenta, however, 
we assume a maximum spin entanglement by choosing the spin part to be one of the Bell states. 
Therefore, we can make a column vector as 




x [D aW ,(W(A 1 , Pl ))* a ,, a ,,(p l ,p 2 )D aW ,(W(A 2 ,p 2 )) 
x D ai CTf ( W (Ai , pi ))* CT »v 2 " (Pi , P2) D a2 a >>> ( W(A 2 , p 2 ) ) 



(19) 



*(Pi,P 2 ) = /(Pi)/(P 2 )Xi 



(20) 
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where /(p) is a normalized function and Xi is is one °f the Bell states 



Xi 



V2 



X3 



V2 



M 



V 1 / 
/o\ 
1 
1 

V / 



X'2 



V2 



V2 



( 1 \ 




V- 1 / 

( \ 
1 

-1 

, , 



(21) 



Substituting the components of (|20p in (|19p . we obtain the final reduced density matrix components 
as 



SI 



\ J ^Pl|/(Pl)| 2 / d 3 P 2 l/(P 2 )| 2 

x [ J D (Tit (e 1 ) J D 4t (e 2 ) J D; it (G 1 ) J D; 2t (e 2 ) 

±D a[t (@ 1 )D a , 2t (G 2 )D* aii (Q 1 )D; 2i (G 2 ) 
±D a[i (e 1 )D a , 2i (Q 2 )D* lt (@ 1 )D^(Q 2 ) 
+D a , i (@ 1 )D^(e 2 )D* aii (@ 1 )D* 2i (@ 2 )}, 



where the upper (lower) sign corresponds to the choice xi (x 2 ) ; an d 



(22) 



(23) 



\ J d 3 Pl |/( Pl )| 2 | d 3 P2 |/(p 2 )| 2 

x [ D CTit (e 1 )D (7 , i (e 2 )^ lt (e 1 )L'^(G 2 ) 
±D a , ii {® l )D a ,j(e 2 )Di i ^e l )D* a2i (e 2 ) 

±D a[t (e 1 )D a2i (Q 2 )D* aii (Q 1 )D; 2t (@ 2 ) 
+D <i {Q l )D^(e 2 )Dl ii {Q l )Dl 2 ^Q 2 ) 



where the upper (lower) sign corresponds to X3 (Xi)- Here the components of -D(Gi) or -D(© 2 ) are 
given by (fTH)) . 

We have suggested that the centroid moves along a circle with a definite momentum, however, 
from a quantum mechanical point of view we imagine that the momentum represented wave packet 
is distributed properly around the momentum of the centroid. In what follows, we choose a suitable 
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form for the distribution function as 



mc 



(p 3 _ g3) 2 



2f3 2 m 2 c 2 



(24) 



which describes a sharp distribution about p 1 = and p 2 = 0, and a Gaussian distribution about 
p 3 = q s with a width specified by /3. 

Let us apply ((24l) in ([22]) and ([23]) and define 



C 



1 



dp exp 



and 



1 f 00 
5 = - — r / dp exp 



(g - qf 
p 2 



(p - qf 

P 2 



cos 0, 



(25) 



sin ©, 



(26) 



3 d 

where p = and q = are dimensionless parameters and is given by 



mc 
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-ere 



-B(r) 



A'frO 



(27) 



:qVq 2 + i\ Vq 2 + i 



vv + 1 + 1 / ' 



according to (|16j) . Then, after doing some manipulations we obtain the final reduced density 
matrices as 



,(/) 



/ 1 ± C 2 ± S 2 

1 T C 2 t5 2 



1 ± C 2 ± S 2 \ 

-1±C 2 ±S 2 



\1± C 2 ± S 



1±C 2 ±S 2 




1 t C 2 t 5 2 







(28) 



i ± c 2 ± 5 2 y 



where the upper (lower) sign corresponds to \i (X4)> an d 



,(/) = - 



/ 1±C 2 t5 2 ^2CS t2C5 -1 t C 2 ± S 2 \ 

=F2C5 It C 2 ± S 2 1=f C 2 ± S 2 ±2CS 

=F2C5 It C 2 ± 5 2 1=f C 2 ± S 2 

\-l^C 2 ±S 2 ±2CS ±2CS 



±2CS 

72 — i — c2 



(29) 



1 ± C T <5 / 



where the upper (lower) sign is used for %2 (X3)- 
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It is interesting to consider a situation in which the centroid falls along a radial geodesic with 
a velocity v, as measured in a local inertial frame. Then we can write v = ce 1 r dr/e° t dt = 
e~ A ^ e B ^ {dr / dt) which leads to q° = ^mc and q 1 = ^mv as the non-zero components of the 
four-momentum of the centroid. Since this is a geodesic motion, the acceleration part of @j van- 
ishes and the curvature part gives only A°i = \ l $ = —cyA'(r)e~ B ( r ) which consists of a boost along 
the 1-axis. In the present case, all the components of w t j are zero as (JBJ) implies. Thus, the Wigner 
rotation (J7J) is trivial and its two-dimensional representation simply becomes D a i a = S a ' a , which 
identically leads to the result q^i = q(v independent of spin entanglements. This means, in the case 
of radial motion the reduced density matrix is invariant and there is no spin decoherence. Of course, 
there exists still a global Lorentz transformation obtained from ([5]) which transforms the initial 
wave packet ([2]) into a boosted frame along the 1-axis via the operator U (Ax(r/, n))® {7(A2(r^, fi))- 

A possible application of the above result is to suggest a general relativistic invariant protocol 
for quantum communication. For an observer falling on a radial geodesic the information stored 
in spins of two entangled spin-^ particles remains perfect. 



IV. SPIN ENTANGLEMENT 



The entanglement between spins as viewed by an observer located at a local frame at the 
final point is obtained by calculating the Wootters' concurrence [3], denoted as C(g). It is 
defined as C(g) = max{0, Ai — A2 — A3 — A4} where AjS are the square root of the eigenvalues, in 
decreasing order, of non-hermitian matrix gg where g = (a y (g> a y )g*(a y (g> cr y ). Note that each Aj 
is a non-negative real number. Then, the entanglement E(g) can be calculated as 

E(,)^(i±4^ (30) 

where 



h(x) = —x log 2 x — (1 — x) log 2 (l — x). (31) 

After doing some manipulations it turns out that the eigenvalues of g^> q"> for all of the matrices 
221) and (i23l are the same and are obtained as 



{i(l + C 2 +5 2 ) 2 ,i(l-C 2 -5 2 ) 2 ,0,o}. 
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Then the concurrence in this case becomes C = C 2 + S 2 . This means that the entanglement for all 
of the states can be calculated by a unique relation as 



W^ lo ( 32 ) 



V. DISCUSSION 

For evaluating these entanglements we need to specify the back ground gravitational field by 
determining the unknown functions A(r) and B{r) in (|27p . The properties of these functions 
affect directly the behavior of the Wig ner rotation in terms of r. In the following we illustrate the 
argument for a charged black hole 



15j |. As we will see below the particular forms of A(r) and B{r) 
for this spacetime provides interesting consequences for the behavior of the entanglement in terms 
of r. The existence of two parameters corresponding to the mass and the charge of this black hole, 
allows us to furnish a detailed argument on the gravitational spin decoherence. This seems to be 
enough to motivate us to consider such a gravitational field. Therefore, we choose 



r 2 

r r 2 



e 2A(r) =e -2B(r) =1 _^ + S_ ^ 



where r s and e relate to the mass and the charge of black hole, respectively. If r s 2 < Ae 2 the 
expression (|33|) will have no real roots. Hence, it follows that the metric Q is singular only at the 
origin r = 0. This singularity is intrinsic and in the present case is called naked singularity, since 
there is no event horizon. The more interesting case occurs when r s 2 > 4e 2 , such that the metric has 
two (or one) other singularities where (|33p vanishes, namely, at r± = (r s ± V r s 2 — 4e 2 )/2. These 
roots are called event horizons. At any event horizon the tetrad (jlip diverges and consequently the 
local Lorentz transformation (j3]) diverges. Therefore, the Wigner rotation angle in ()27p must 
diverge at any horizon. 

Substituting ([33]) in (f2T|) . the angle can be written as 




2z 2 y/z 2 -z + g 



V Vp 2 + 1 + 1 / 
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FIG. 1: Plot of E for fixed £, z, r and /3 but varying q. As q increases more spin dccohcrence occurs and so 
the entanglement descends. 

where t s is a constant proper time during which a photon rotates once on a circle of radius r s , and 
£ 2 = -^7 and z = — are dimensionless parameters. 

For £ 2 < 1/4, we encounter two event horizons located at z± = (1 ± yl — 4£ 2 )/2. As (|34D 
implies, in the region between the horizons, that is z_ < z < z+, the angle B becomes imaginary 
and our argument fails there. So we consider only the region z > z+. (The region z < z_ is not 
physically interesting in our argument). In this case has two zeros however only one of them 
Z\ = (3 + — 32£ 2 )/4 is accessible. For £ 2 = 1/4, only one event horizon occurs which is located 
at z = 1/2 and has an accessible zero at z = 1. In the case of naked singularity £ 2 > 1/4 and 
there is no event horizon. Then is real everywhere and it is singular only at z = 0. There exist 
two accessible zeros for at zi 2 = (3 ± v / 9-32£ 2 )/4 provided £ 2 < 9/32, other wise it has no 
zeros. 

Substitute given by (|34|) in ([25]) and (|26|) , it turns out that after integrating, the entanglement 
E given by (132[) depends formally to the parameters £, r, /3, 2 and There is no analytical solution 
for the integrals and so we attend to a numerical approach for solving them. In the following we 
indicate the behavior of E in some graphs all obtained numerically. 

In Fig. [H the entanglement is sketched versus a q for £ 2 = 0.265, r = 5r , z = 1.6 and (3 = 1. 
We see that by increasing the angular velocity of the centroid on the circle, E descends uniformly 
from 1 to an asymptotic value of zero. Such a behavior is justified because, for a constant r, by 
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FIG. 2: The same plot of figure I but for a larger value of j3, that is a wider momentum distribution. As q 
increases the curve again descends but by performing aperiodic oscillations. 

increasing q, the centroid travels more on its circular trajectory in the gravitational field and so 
more spin decoherence occurs. It is remarkable that increasing the width /3, has a meaningful effect 
on the behavior of E in terms of q. This point is shown in Fig. [2] which is plotted for /3 = 4 and 
the other parameters as those of Fig. [TJ The entanglement again diminishes but by performing an 
aperiodic oscillation. This can be justified by noting that as (3 increases, the exponentials in the 
integrals (|25p and (|26p approach the unity , so the cosine or sine terms play a more dominant role. 

One may also plot E in terms of t/tq where r is the proper time during which the centroid 
moves in the gravitational field. Fig. [3] shows such a graph. Again as is expected, increasing r 
produces more spin decoherence and the entanglement falls. 

More interesting here is to plot the entanglement in terms of the radius of the circular orbit 
characterized by z. Fig. [J] shows E versus z for £ 2 = 0.16, $ = 1, q = 0.6 and t = 5r . In this 
case there exist two horizons at z = 0.2 and z = 0.8, however the lower limit of z is taken to be 
the outer horizon z = 0.8. At any horizon © diverges, then the cosine and sine terms in (|25p . (|26p 
oscillate rapidly yielding a zero value for the integrals. Then the entanglement (j32|) is completely 
destroyed, leading to a fatal error in quantum communication at the horizon. On the other hand, 
the accessible zero of O occurs at z = 1.24 where the Wigner rotation (I17p becomes the identity 
operator. On such a circle the curvature part and the acceleration part of @ cancel out each other, 
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FIG. 3: Plot of E for fixed £, z, q and j3 but varying r. Increasing r also causes to occur more spin 
decoherence and so the entanglement diminishes. 

leading to a zero value for 0. Thus, for the corresponding circular orbit E remains the unity, that 
is, there is no decoherence for the spin Bell states. This result is independent of values of q, j3 and 
r. As z grows, the curve takes a minimum at z = 2.25 which determines a circular orbit on which 
the effects of acceleration and curvature add to produce a maximum spin decoherence. Note that 
at z — > oo, both the curvature part and the acceleration part of vanish. Hence, as (|34p shows, 
lim^oo = and the curve asymptotically approaches to 1. 

In the case of naked singularity only diverges at z = 0, where, the entanglement vanishes. 
Fig. [5] and Fig. [6] show E in terms of z for this case. However, the curve in Fig. [5] is depicted 
for £ 2 = 0.265 so has two accessible zeros at z = 0.57 and z = 0.93 where the Wigner rotation 
becomes the identity operator and so the entanglement remains intact. There is no spin decoherence 
on these two circles, independent of the values of q, (3 and r. Also the curve takes two minima at 
z = 0.65 and z = 1.96 indicating two circles of maximum decoherence. The curve asymptotically 
turns to the unity. In Fig. [6] the curve is sketched for £ 2 = 0.5 and then there is no zeros for 0. By 
increasing z, the entanglement grows from zero at z = and reaches asymptotically to the unity. 

That on some circles E remains the unity may suggest a possibility for a perfect quantum 
communication around a charged black hole. The properly entangled spins of the system can 
preserve the stored information during the motion on such circles. 
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7* 1 1.5 2 2.5 3 3.5 4 

outer ' ^ 
horizon 

FIG. 4: Plot of E for fixed £. r, g and /3 but varying z. Here £ is chosen such that there exists two event 
horizons for the spacetimc. z begins from the outer horizon where E vanishes. There is a peak which specifics 
the radius of a circular orbit on which no spin decoherence occurs and the entanglement remains the unity. 
Also, the curve takes a minimum which describes a circular orbit with maximum spin decoherence. The 
curve asymptotically reaches to the unity. 



0.3- 
0.2 : 
0.1 - 




FIG. 5: The same plot of figure 4 but for the case of naked singularity. In the vicinity of singularity E 
vanishes. £ is chosen such that has two zeros at which the peak of curve occurs. The peaks designate two 
circles on which the entanglement remains intact. Also the curve takes two minima indicating two circles of 
maximum decoherence. The curve asymptotically turns to the unity. 
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FIG. 6: Plot of E versus q for the case of naked singularity, but £ is chosen such that no zeros occur for 0. 
Then the curve grows from zero at z = and reaches asymptotically to the unity. 



VI. FRAME TRANSFORMATION 



Consider a local transformation between two instantaneously coincident inertial frames, as 

ea" — > e/ = T a b (x)e b », (35) 

where T a c {x)Tb d {x)r] C( i = rj a b. Notice that the definition (pQ) remains intact under this tetrad 
transformation which can arise from a coordinate transformation x^ — > x^. However, as (j3|) 
implies, the new inertial observer who uses ej 1 obtains different values for the components X a b(x), 
and regarding ([8]), we recognize that the Wigner rotation is not invariant under the transformation 
(|35p . Consequently, the spin entanglement is not frame independent. To illustrate the situation, 
let us simply consider the Schwarzschild black hole by letting e A ^ = e~ B ^ = */l — ^ in the 
metric Q. Consequently, the tetrad (fTTj) reduces to 

eo l = yL= , ei r = ^=L=, e 2 6 = -, e ^ = — (36) 
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which is used by a static local observer. Note that this is singular at the horizon r = r s . Then, 
the local Wigner rotation elements (|16p reduce to 



w 1 3 



-w 3 i 



1 _ 

1 2r 



qp 



vV+T + i 



(37) 



which can be used for finding the elements of the Wigner rotation operator and then calculating 
the entanglement (|32j) . These elements diverge at the horizon and so the entanglement will be 
completely destroyed there, leading to a fatal error in quantum communication. Also they vanish 
at the circle r = ^ where the entanglement will remain intact. Now, let us transform the present 
spherical coordinates (t, r, 6, cf>) to, for instance, the Kruskal coordinates (T, R, 8, <f>) via the relations 



R 2 - c 2 T 2 



cT / ct 

—— = tanh 

R V 2r, 



(38) 



In the Kruskal coordinates the Schwarzschild metric becomes 



ds 2 = ^e~^ {-c 2 dT 2 + dR 2 ) + r 2 (d8 2 + sin 2 ^# 2 ) 



(39) 



where the radial coordinate r is now interpreted as a function of T and R. Accordingly, we choose 
a new tetrad Ca as 



S T 



1 IT 



c V r 



e 2r ", e 1 



R 



r — 1 

e 2r s e = _ e 



1 



r sin ( 



(40) 



which is not singular at the horizon. According to (j35H . we find that the this tetrad is related to 
the static tetrad (|36h by a local transformation 



T a C {x) 



R 



2 V r s (r-r s ) 2 V r s (r-r s ) 



c / e r s 



TOO 



i 



2 y r s (r-r s ) 2 y r 3 (r-r s ) 



V 



o 
o 








R 

1 
1 



(41) 



This transformation resembles a Lorentz boost along the 1-axis and so the new local inertial 
observer falls into the black hole when T > 0. Now, repeating the steps of Section 2 with the 
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tetrad (flOj) . we obtain, instead of ([37|) . the local Wigner rotation elements 

^-^--i-M{ l+ y)^{^-7whr)' (42) 

where q and p relate to the four-momentum of the centroid and the four- momentum of particles, 
respectively, that are chosen by the falling observer. Since these elements are not singular at the 
horizon, the spin entanglement can have a nonzero value and falling observers just on the horizon 
can perform a quantum information task with a finite error. Furthermore, they can do this beyond 
the horizon r < r s , until the physical singularity r = 0. Therefore, the features of quantum 
communication in a curved spacetime using the entangled spins depend on the tetrad that is used 
by the observer. 

VII. CONCLUSION 

In this paper we discussed the effect of a gravitational field on the spin entanglement of a 
system of two spin-| particles moving in the gravitational field. We described the system by a two- 
particle Gaussian wave packet represented in the momentum space with a definite width around the 
momentum of a mean centroid. We considered the mean centroid for describing the motion of the 
system in the gravitational field which was a classic field in our approach. Both the acceleration of 
the centroid and the curvature of the gravitational field caused to produce a Wigner rotation acting 
on the wave packet. We confined the argument to spherically symmetric and static gravitational 
field such as a charged black hole. At a given initial point in the curved spacetime, the momenta 
of the particles were taken to be separable, while the spin part was chosen to be entangled as the 
Bell states. By calculating the reduced density operator, we focused on the spin entanglement of 
the system. 

For any spherically symmetric and static spacetime, when the system, no matter of spin and 
momentum entanglements, falls along a radial geodesic, the Wigner rotation is trivial and so the 
spin density matrix remains intact. Moreover, for a charged black hole, depending on the mass 
and the charge of black hole, there can be circular paths with determined radii on which the spin 
entanglement is invariant. For instance, it always remains the unity, provided spins be entangled 
initially as one of the Bell states. This result is a consequence of cancellation of the acceleration and 
the curvature on such circles. The spin entanglement descends to zero by increasing the angular 
velocity of the centroid as well increasing the proper time during which the centroid moves on its 
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circular path around the center of charged black hole. 

It is remarkable that the features of quantum communication in a curved spacetime using the 
entangled spins depend on the tetrad that is used by the observer. 

As an extension of this work, one can provide the initial state with a specified spin entanglement 
and also a specified momentum entanglement, then discuss the final spin entanglement, as well as 
the final momentum entanglement. 
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